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Abstract 

Transformation models provide a popular tool for regression analysis of censored 
failure time data. The most common approach towards parameter estimation in these 
models is based on nonparametric profile hkelihood method. Several authors proposed 
also ad hoc M-estimators of the Euclidean component of the model. These estimators 
are usually simpler to implement and many of them have good practical performance. 
In this paper we consider the form of the information bound for estimation of the 
Euclidean parameter of the model and propose a modification of the inefficient M-- 
estimators to one-step maximum likelihood estimates. 



1 



February 2, 2008 



2 



1 Introduction 

The proportional hazard model, originating in Lehmann [17] and Cox [8], provides the 
most common tool for regression analysis of failure time data. The simplest version of 
this model assumes that a failure time T and a covariate Z satisfy the interrelationship 

h{T)^-P'^Z + e (1.1) 

where h is an unknown continuous increasing function mapping the support of T onto 
the real line, /3 is an unknown regression coefficient and e is an error term independent of 
the covariate Z and having extreme value distribution with density f(x) = exp[— e^]. 
Alternatively, if fi is the marginal distribution of the covariate and A{t, z) the cumula- 
tive hazard function of the conditional distribution of the failure time T given Z = z, 
then the model stipulates that 

'A{t,z) = T{t)e^^' l^a.e.z, (1.2) 

where j3 is an unknown regression coefficient and F is an unknown continuous increasing 
function, F(0) = 0, mapping the support of the failure time T onto the positive half 
line. 

During the past two decades several authors proposed generalizations of (1.1)-(1.2) 
to semiparametric transformation models specifying the interrelationship between the 
conditional hazard function A{t, z) and the transformation T{t) as 

A{t,z) ^ A{V{t),e\z) na.e.z, (1.3) 

where {A{-,9\z) : 9 e 6} is a family of conditional cumulative hazards dependent on 
a finite dimensional parameter 9. Following Bickel et al. [4], this family is referred to 
as the "core model". Common choices of (1.3) include scale regression models with 
core model derived from distributions with decreasing hazard rates, such as the frailty 
distributions with finite mean. In particular, the proportional odds ratio model has 
gained much popularity as a competitor to the proportional hazard model [1, 9, 10, 15, 
18, 21, 23]. Core models with increasing or non- monotone hazards were considered in 
[5, 6, 7, 9, 12]. 

In this paper we consider estimation of the parameter 9 based on an iid sample of 
right-censored failure times. For purposes of analysis of the odds ratio model. Murphy 
et al. [18] proposed to use nonparametric profile likelihood method. The approach 
taken was similar to the classical proportional hazard model. The model (1.3) was 
extended to include all monotone functions. With fixed parameter 9, an approximate 
likelihood function for the pair {9, F) was maximized with respect to F to obtain an 



February 2, 2008 



3 



estimate Tns of the unknown transformation. The estimate Tno was shown to form 
a step function placing mass at each uncensored observation, and the parameter 6 
was estimated by maximizing the resulting profile likelihood. Under certain regularity 
conditions on the censoring distribution, the authors showed that the estimates are 
consistent, asymptotically Gaussian at rate y/n. They also proposed to estimate the 
standard errors of the regression coefficients by numerically twice-differentiating the 
log-profile likelihood. The approach was generalized to other transformation models 
whose core models have decreasing hazards in [16, 21, 22]. 

In this paper we take a different approach towards construction of efficient estima- 
tors of the parameter 6. To motivate it, let us recall [4] that in a regular parametric 
model V = {Pe,r] : 9 E Q,ri E H}, the asymptotic variance of any regular estimator of 
the parameter 9 satisfies the bound 

yar[V^{9-9)] > [h,{9,ri) - Ii2{9,v)l22{9,vr'l2i{9,v)]-\ 

where 

r(n (In{9,v),h2{9,v) 

'^^'^^ [h,{9,7l),h2{9,7j) 

is the Fisher information matrix with entries 

In{9,v) = Eie{X,9,rif^ l22{9,v) ^ Ee,{X,9,rif^ 
1,2(9, v) = EieiX,9,r))i,{X,9,r)f ^ I2ii9,r)f. 

Here £$ and £^ represent score functions corresponding to the two parameters. Alter- 
natively, 

var[v^(^-e)] > [Et{X,9,rj)e*{X,9,rif]' 

where £* is the efficient score function for estimation of 9. The function £* is the 

(componentwise) projection of the vector of scores ig onto the orthocomplement of 
the nuisance tangent space Vrj spanned by all scores of the nuisance parameter rj. To 
estimate the parameter 9, we may consider solving the score equation 

n ^-^ 

1=1 

where ri{9) is an estimate of rj obtained "for each fixed 9" in the parameter set 0. To 
be more precise, we assume that V is a submodel of a larger family of distributions 
Q, and there exists a parameter rj : Q x ^ 7Y, with the property rj{Q, 9) = t], 
whenever Q = Pg jj G V. We require the estimate r]{9) to be asymptotically unbiased 
for estimation of r]{Q, 9) in the larger model Q. If the equation (1.4) has a consistent 
solution 9, then 9 is an efficient estimate of the parameter 9 under additional conditions 
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on r]{9)) given in [4, Ch.7.7]. These amount to the assumption that i*{-,9,r]{9)) is 
a consistent estimator of the efficient score function and the bias of this estimator 
converges in probabihty to at a rate faster than ^/n. 

Turning to transformation models, in the case of uncensored data Bickel [2] , Bickel 
and Ritov [3] and Klaassen [15] used invariance of the model with respect to the group 
of increasing transformations to show that the efficient score function for estimation 
of 9 is given by a nonlinear rank statistics. Its form was derived using Sturm-Liouville 
theory. In the case of censored data, several authors verified existence of y/n estimators 
of the unknown transformation and used them to construct ad hoc estimators of the 
parameter 9. Whereas these estimators are inefficient, many of them have good prac- 
tical performance [5, 6, 7, 23] and are simpler to implement than the proffie likelihood 
method. They also apply to a wider class of transformation models. 

In section 2, we assume the so-called " non- informative censoring model" and con- 
sider core models whose hazard rates are supported on the whole positive half-line, 
and are finite and positive x — 0. Under a certain integrability condition, we derive 
the form of the information bound and efficient score function for estimation of the 
parameter 9. In section 3 we verify the integrability condition in the special case of the 
generalized odds ratio and the linear hazard regression models. In Section 4 we con- 
struct a class of Z estimators of the parameter 9. To study its properties, we shall make 
the assumption that the censoring distribution has support contained in the support of 
the failure time distribution and its upper point forms an atom. Under mild regularity 
conditions, we show that the proposed Z-estimators have an asymptotic distribution 
not depending on the choice of the estimator of the unknown transformation. As a 
by-product, we also show that the parameter 9 can be efficiently estimated by solving 
an equation of the form (1.1) or by means of one-step MLE. 



2 Information bound 



2.1 Martingale identities 

Throughout the paper we assume that the triple {X, 5, Z) represents a nonnegative 
withdrawal time {X), a binary withdrawal indicator {5 — 1 for failure and S — ioi 
loss-to-follow-up) and a vector of covariates (Z). The triple {X,S,Z) is defined on a 
complete probabihty space (f), ^, P) and X — T AT' and 5 = 1{X = T), where T and 
T' represent failure and censoring times, respectively. We assume that T and T' are 
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conditionally independent given Z. 
We denote by 

'A{t,Z) = I P{X ^du,8^\\X>u,Z), 
Jo 

A{t,Z) = P{X edu,5^0\X >u,Z), 
Jo 

the conditional cumulative hazard functions of the failure and censoring time and set 

N(t) = l{X<t,5^1), Nc(t) ^1(X <t,5 ^0), 

A{t) = f Y{u)A{du,Z), Ac{t)= [ Y{u)Ac{du,Z), 
Jo Jo 

where Y(u) = 1(X > u). Then M(t) = N{t) - A{t) and Mc(i) = N(t) - Ac{t) 
form mean zero martingales with respect to the sclf-exciting filtration {!Ft : ^ > 0}, 
- cr{l{X < t)d, 1{X < t)X, 1{X <t){l- S), Z}. 

Let To = sup{t : EY{t) > 0} and let Q and Qc denote the joint subdistribution 
functions 

Q{t,z) = P{X <t,S = l,Z <z) and Q^t, z) = P{X < t,S = 0, Z < z). 
Then the failure and censoring counting processes satisfy 

fTO 

h{u, Z)M{du) E Ll{P) iff /i e L2(g), 

/■TO 

/ h{u,Z)Mc{du) eLl{P) iff heL2{Qc). 
Jo 

Assuming that the function A is continuous, the two martingale processes are orthog- 
onal. Nan, Edmond and Wellner [19] showed also that any function b e can be 
represented as a sum 

/■TO pro 

b{X, 5,Z)^ / Ri[b]{u, Z)M{du) + / R2[b\{u, Z)M^{du) + E[b\Z\, 
Jo Jo 

where Ri[b] e L2{Q) and R2[b] G I/2(<5c) are given by 

Ri[b](x,z) = b(x,l,z) - E[b(X,5,Z)\X >x,Z ^ z], 
R2[b]{x,z) = b(x,0,z) - E[b(X,5,Z)\X >x,Z ^ z]. 

The three terms in this representation are orthogonal. 



/ 

Jo 
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2.2 Assumptions and notation 

To derive the form of the information bound and efficient score function in the trans- 
formation model, we make the following regularity conditions on the core model. 

Condition 2.1 (i) The parameter set C -R'^ is open, and the parameter 9 is 
identifiable in the core model: ^ 7^ ^' iff A{-, 9, z) ^ A{-, 9', z) for //-a.e. z. 

(ii) For // almost all the function A{-,9,z) has a hazard rate a{-,9,z) supported 
on the whole positive half-line and there exist constants < mi < m2 < 00 such 
that mi < Q!(0, 6*, z) < m2 for ii-Si.e. z and all 6* e 0. 

(iii) The function i{x, 9, z) — logQ;(x, 9, z) is continuously differentiable with respect 
to both X and 9. 

(iv) The family {a{x, 9,z) : 9 & 0} forms a regular parametric model, in particular 
information is finite and positive definite. 

The derivatives of i{x, 9, Z) are denoted by 



Under assumption of transformation model, the true distribution P of (X, 5, Z) 
belongs to the family V — {P{e,r]) : 9 E Q,r] E H} where 9 is the Euchdean parameter 
of interest and rj = (F, Ac, n) is the nuisance parameter corresponding to the unknown 
transformation F, the unknown cumulative hazard function Ac{t; z) of the conditional 
distribution of the censoring time given the covariate Z and the marginal distribution 
II of Z. We make the following regularity conditions. 

Condition 2.2 Let tq = sup{t : P{X >t)>0} and tf = sup{t : P{T >t)>0}. 

(i) The distribution /i is nondegenerate. 

(ii) The parameters A,, and fi are noninformative on {9, F). 

(iii) F is an increasing continuous function, F(0) = and limt^Tp r(^) = 00. 




d 



d 

dx 



i{x,9,Z). 
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(iv) If To < Tp then tq is an atom of the marginal survival function of the censoring 
time. 



The condition (iii) refers to the unobserved model corresponding to uncensored 
data. If To < Tp, then the condition (iv) serves to ensure that the the transformation 
can be estimated consistently only within the range [0,ro]. We do not know whether 
there exists a consistent estimate of the parameter ^ if To, To < Tp, is a continuity point 
of the censoring distribution. 

Using conditions 2.1-2.2 and the assumption of the conditional indpendence of 
failure and censoring times, the score operators for the four parameters are 



[r{u),e,z)M{du), 



TO 



TO 



g{u) + e'{r{u),e,z) 

b{u,Z)dMc{u), 



mix,s,z) = 

e2[g]{X,d,Z) = 

ismx,5,Z) = 
h[c]{X,S,Z)) = c{Z), 
where c{Z) e L'^in), b{u, Z) e L2{Qc) and 

g{x) + e'{r{x),e,z) [ gdvei^iQ) 

Jo 



gdr 



M{du), 



(2.1) 



The tangent spaces for the four parameters are Vi = = 1, .... 4, where [o] denotes 
the closed linear span of the set a in L2(P). The spaces Vi,i = 2,3,4 are mutually 
orthogonal, and so are the spaces Vi,i — 1, and i = 3,4. The nuisance tangent space 
is V^^V2 + Vz + Va. 



In the case of the proportional hazard model, 

Z(i,z) = r(i)e^"^ (2.2) 
Sasieni [20] and Nan, Edmond and Wellner [19] showed that the martingale oparator 

W)= r {f{x,Z)-E[f{X,Z)\X^x,5^l])M{dx), feL^iQ), 
Jo 

satisfies (i) U{f) ± V, in L^iP); (ii) n(/|P^^) = U{R,{f)) for / e Ll{P) and (iii) 

"P^ = {U{f) : f e L2{Q)}. Here Il{-\V^) denotes the projection onto V^, the ortho- 
complement of the nuisance tangent space. These properties entail that the efficient 
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score function for estimation of the parameter 9 in the model (2.2) is given by 
f{x,Z) = Z. 

To obtain an extension of this result to the transformation model, we shall use the 
following notation. Firstly, we shall find it convenient to denote the "true" parameters 
(^,r) as (^o,ro). Further, let 

s[l]{u,T,e) = EY{u)a{T{u),e,Z) 
s[f\{u,V,e) = EY{u)f{u,Z)a{V{u),e,Z), 

e[/]Kr,^) = ^Kr,^) 

and 

cov[/i, h] = e[/i/J] - e[mf2f , var[/] = cov[/, /] 

Note that if failure and censoring times are conditionally independent given Z, then 
under the assumption of the transformation model, the conditional distribution of Z 
given X — t and 5 — 1 has form 

P{ZeB\X = t,6=l) = I I G,{t,z)F{ro{t),eo,z)a{ro{t),eo,z)fi{dz) 

= e[/](i,ro,^o), 

where f{t,Z) = 1{Z G B), F{x,9o,z) = [exp —A{x,9o, z)] is the survival function of 
the core model and Gc is the conditional survival function of the censoring time. More 
generally, 

e[f]{t,ro,eo) = E[f{X,Z)\X = t,6 = 1] 

and similarly, var[/](t, Fq, 6^0) and covo[/i, /2](t, Fq, ^0) are conditional variance and 
covariance operators. From section 2.1, we also have 



fiX,Z)eL2{Q) iff E jY{u)fiu,Z)airoiu),eo,Z)roidu) 
Using domintated convergence theorem, 

fix,z)eL2iQ) iff s[f]iu,ro,9o)eL^{ro). 

By noting that 

f ENjdu) 

the square integrability condition reduces to 

f{X,Z)eL2{Q) iff e[f]{u,ro,eo)eL,{Q). 



< 00. 
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Next, with some abuse of notation, we shall write s[£']{u,r,9),e[£']{u,r,9) and 
var[f](M, r, 6) whenever f{X, Z) = f(r(X), 6, Z). Let 

C{t) = [ {s[l]{u,r,e))-^ EN{du), 
Jo 

B{t) = [ vaT[i']{u,T,e)EN{du), 
Jo 

V{u,t) = exp- /" s[e']{u,r,0)C{du), (2.3) 

J u 

ptAt' 

K{t,t') = / C{du)V{u,t)V{u,t') . 
Jo 

Finally, let D[f],D[f]{t) = D[/](i,r,^), denote the solution to the linear Volterra 
equation 

D[f]it) = - fs[f]{u,T,e)C{du)- f D[f]{u-)s[e']{u,T,e)C{du) 
Jo Jo 

The equation has a unique locally bounded solution given by 

D[fm = - f s[f]{u,T,e)C{du)V{u,t). 
Jo 

We shall need the following integrabihty condition. 
Condition 2.3 Let 

k{to)^ [ [ C{du)V{u,tfB{du) 

J J0<U<t<T0 

and suppose that ac(to) < oo. 



Lemma 2.1 If the condition 2.3 holds then K e L2{B B) and D[f] e L2{B), for 
any function f{X,Z) such that e[/^] e Li{Q) 



This lemma can be verified using Cauchy-Schwartz inequality and dominated con- 
vergence theorem. If failure and censoring times are conditionally independent and 
the transformation model is satisfied by {9,T) = (6'o,ro), then the second part of the 
lemma holds for any f{X, Z) e L2{Q). 
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2.3 Information bound 



Unless this leads to confusion, we suppress the dependence of the functions s[/](-, T, 6), 
e[/](-,r,^) and the corresponding variance and covariance operators on (T,^). We 
write So[/], eo[/], varo[/] and covo[/i,/2] whenever the failure and censoring times are 
conditionally independent given Z and the transformation model is satisfied by (F, 6) — 
(Fo, ^o)- The functions Cq,Bq,Kq, Vo and Dq are defined analogously. 

Further, for any ip e L2{B) and f{X, Z) such that e[/^] e Li{Q), set 

p[f,cp]{u)^coy[f,e']{u)-cp{u)y^r[e']{u) (2.4) 

and let 

Wf(t, e, z) = fit, z) - e[f]{t) - (/(F(t), e, Z) - e[e]{t)) v*{t) 

^l^^l nt,s)p[f,^*]{s)EN{ds), (2.5) 

where (p* is the solution to the Fredholm equation 

(p*{t) ^ -D[f]{t) - K{t,u)(p*{u)B{du)+ K{t,u)cov[fJ']{u)EN{du). (2.6) 
Jo Jo 

If var[f ] = 0, then Wf{t, 9, Z) = f{t, Z) - e[f]{t). In addition, setting ^p* = ^* + 
the equation (2.6) simplifies to 

/"TO 

V'(0 - A / K{t,u)ip{u)B{du) = r]{t), (2.7) 
Jo 

where A = — 1 and 

rj{t) = J K{t,u)p[f,-D[f]]{u)EN{du), 
p[f,-D[mu) = coy[f,e'](u)+Y8.r[£']{u)D[f]{u) 

If var[£'] ^ then the kernel K is symmetric, positive definite and square integrable 
with respect to B. Therefore it can have only positive eigenvalues. For A = — 1, the 
equation (2.7) has a unique solution given by 

r{t)^v{u)- A{t,u,-l)r){u)B{du), (2.8) 
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where A{t, u, A) is the resolvent corresponding to the kernel K. The resolvent equations 

K(t,u) = A{t,u,X)-X A{t,w,X)B{dw)K{w,u) 

Jo 

= A{t,u,X)-X K{t,w)B{dw)A{w,u,X), 
Jo 

applied with A = — 1 imply 

r{t)^^*(t) + D[f]{t) = r A{t,u,-l)p[f,-D[mu)EN{du). (2.9) 

^0 

If var[£'] ^ but p[f,—D[f]] = 0, then the solution to this equation is ip* = 0, or 
equivalently, (p* = —D[f]. 



Define 



u{f,e)^ rW{f,0,Z){u)M{du) 
Jo 



Proposition 2.1 Suppose that the conditions 2.1-2.3 are satisfied with {6, T) = {6q, Fq). 
Then (i) U{f,eo) ± in L^iP); (ii) Il{f\V^) = U{R,{f),9o) for / e L0(P) and (in) 
V^ = {U{f,9o):feL,{Q)}. 



Proof . Set 



1 



To 



so[l]{t)jt 

Then ^ 

Jo 

and 



Vo{t,s)po[f,<p*]{s)EN{ds). 



PTO pro 

U{fM= / f{t,Z)M{du)-t2W\^ / [f{t,Z)-g*{t)-e{V^{u),eo,Z)ip*{u)\M{du). 
Jo Jo 

Using Cauchy-Schwartz inequality and the condition 2.3, it is easy to verify that ^2[fi'*] ^ 
L^{P) for any function /(X, Z) G L2{Q). Therefore part (i) of the proposition will be 
verified if we show that E{U{f, 6*0)^2(5']) = for any function g satisfying th condition 
(2.1). 
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Put 4^ = h{g) + h{g), where 



h{g) 



g{u)dM{u), 



h{g) - r kroH,^o,^) / gdv, 

Jo L Jo 



M{du). 



Then 

EU{fMh{g) 
EU{fMh{g) 



/ git) 

Jo 



TO 



Vo{t,s)po[f,ip*]{s)EN{ds) 



Toit), 



TO 



Po[f,^*m / gdVo 



dEN{t) 



TO 



eoK](t) / Vo{t,s)po[f,ip*]{s)EN{ds) 



gdVo ] dTo{t). 



Using 



'Po{t,s) 



so[i']dCo{du)Vo{u,s) = l- / eo[i']{u)dTo{u)Vo{u,s) 



ior t < s, and applying Fubini theorem, we have 



u 



EU{fMh{g) 

TO 

g{u 

TO 

g{u 

TO 

g{u 

To 



giu) 



TO 



Po[f,<p*]{t)dEN{t)j dToiu) 



TO 



eo[^'] W i^J^ Vo{t, s)po[/, ip*]{s)EN{ds) ) dVoit) ) dVoiu 
Po[f,^*mdEN{t)^ dVoiu) 



To 



eo[i']{t)dTo{t)Vo{t,s) ] po[f,^*]{s)EN{ds) ] dTo{u) 

u \J t 

" Po[f,'P*]it)dEN{t))dro{u) 



[Po{u,s) - l]po[f,ip*]{s)EN{ds) dVoiu) = -EU{f,9o)Ii{g). 



pTQ 

= / 9( 

Jo 

/•To 

+ / g{u 

Jo 

This completes the proof of part (i). Part (ii) and part (iii) follows in the same way as 
in Nan, Edmond and Wellner [19]. □ 

Since U{f, ^o) is a martingale operator, we have 
Eo(/,^o) = EU{f,9of^ 
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= E / Wj{u,eo,Z)Y{u)a{V^{u),e^,Z)V^{du). 



The matrix Eo(/, ^o) satisfies Eo(/, ^o) = ^o) + Sal/, ^o), where 

/•TO 

Si(/,^o) = / VB.Y^[f-e^*]{u)EN{du), 
Jo 



Co{du) 



ro{u,t)po[f,^*mEN{du) 



The conditional variance function is identically equal to 0, if and only if 

fit, Z) = £'{ro{t), 00, Z) [ hdVo + a{t) e L2{Q) 

Jo 

and a = h. For any function of the form (2.11), Fubini theorem yields 

Do[f]{t) = - / a{u)droiu)Voiu,t) 
Jo 

~ lo Uo ^'^ro)^o[^'K«)'^ro(M)Po(«,t) 

= —w{t) — / hdVo, 
Jo 

w{t) = / [a- h]{u)dro{u)Vo{%t). 
Jo 



Hence (2.6) reduces to 



[g*-h]dTo+ / Mt,u)[ 



TO 



{g* - h)dro 



Bo^du) = w{t). 



(2.10) 



(2.11) 



If a = h then the right-hand side of this equation is identically equal to 0, and corre- 
spondingly, g* — h. Otherwise, we obtain 

^p*{t)= I hdTQ + w{t)- I A{t,u,-l)w{u)Bo{du) 
Jo Jo 

and 

/•TO 

varo[/-£V*](i)=varo[f](i)[«;(i)- / A{t,u, -l)w{u)Bo{du)] 

Jo 

If varo[^'] ^ 0, then the right-hand side is identically equal to if and only if 

/•TO 

w{t)^ A{t,u,-l)w{u)B{du). 
Jo 

In this case, the resolvent equations imply that w = so that a = h. 
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3 Examples 



In this section we verify the square integrabihty condition in two models. In particular, 
we show that the information bound applied to both censored and uncensored data. 



Example 3.1 Generahzed odds ratio model. The survival function of the core model 
is given by 

F{x,e,Z) = [l+'qe^^^x]-^-^''' for 77 > 0, 
= exp[— e^^^x] for 77 = 0. 

The proportional hazard model corresponds to the choice 77 = and the proportional 
odds model to 77 = 1. The hazard function of the core model given by 

a{x,e,Z) = e'^^^[l + T]e'^^x\-\ 

We have 

e{x,e,Z) = -^e^^^[l + 77e^^^x]-^ 

If |Z| < do and < e^^^ < ^2, then 

\t{x, e,Z)\ < do[l + rjdix]-^ < do 
r)di[l + rjdix]-^ < -£'{x,e, Z) < r]d2[l + r)d2x]~^ 

Note that —i'{x, 9, Z) is an increasing function of rje^^^ . 

Suppose that (^o,ro) is the "true" parameter of the transformation model. The 
preceding bounds imply 

Vd2 ,^ , . l + 7?o?2ro(t) 



Vo{u,t) < exp [ ^j^^ . dVoiv) 



l + rid2To{u)' 



Next recall that if 6'^ is a random variable with finite variance and distribution function 
H, then 

Var(t/) = i y Jiy^i- W2fH{dwi)H{dw2) (3.1) 
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By noting that 



[e'{ro{t),eo,z,)-e'{ro{t),eo,z2)f^ 

[1 + r;e^^iro(t)]2[l + r7e^<T-^ro(t)]2 
< rf[d2-d^f[l + r,d^Vo{t)]-\ 

and using (3.1), we obtain var[£'](t, Fq, ^o) = 0(l)[l+7yciiFo(i)]~^. Since [l+7yci2Fo(t)] < 
((i2/'ii)[l + ^?'iiFo(t)], we see that 



k{t,) = 0(1) j 



dTr)(v) 



s[i]{v,ro,eo) [i + r)d2ro{v)y^ 



s[l]{t,To,eo)dTo{t) 
{l + r|d^^o{t)y 



In the generahzed odds ratio model, s[l](t, Fq, ^o) is a decreasing function of t. There- 
fore 



s[1](m,Fo,^o) 



< 1 



and 



k{to) < k{oo) = 0(1) + rid2w))-\l + r]diw)-^dw = 0(1) 



The bound is vahd for any distribution of the censoring times. We obtain a better 
bound on the constant k(oo) in the case of the so-called Koziol-Green censoring model. 
The conditional survival function of the censoring time given the covariate is of the 
form 

Gc{t\z)^F{ro{t),9o,zy, a>0. 
The choice of a = corresponds to the case of uncensored data. We have 



so that 



d^[l + r;(i2Fo(t)]-[^+(^+")/''l < s[l](t, Fo, ^o) < d2[l + vdiTo{t)]-^'+^'^'^ 



)/v] 



Hence 



^(oo) = 0(1) 



(l + r;4w)(^+''^/''-Mw 



;i + 77(iix)-t3+(i+«)/'?](ix 



0(1) / {l + r]dix)-^dx = 0{1). 
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Example 3.2 Linear hazard model has failure rate 

a{x, 9\z) = ae{z) + xh0{z). 

We assume that mi < ao{z) < m2, m'l < he^x) < m'2 for some finite positive constants 
mqjTn'qjq = 1,2. In addition, the functions ag{z) and have bounded derivatives 
with respect to 9. We have 

I'ixAz) ^'^'^ 



ag(z) + Xbg(z) 



Set di — m[/m2 and ^2 = m'2/mi and suppose that (^O) Tq) is the true parameter 
of the transformation model. Using a similar algebra as in the example 3.1, we can 
show that 

Vo{u,t) < \ '"^'^f'^^} var[f](t,ro,^o) <0(l)[l + (iiro(t)]-l 



Hence 

«:(ro) = 0(1) 
For V < t, we have 



1 + diToit) 

* dVoiv) 



TO 



s[i]{v,ro,eo) 



[l + diTo{v)f 



s[i]{t,ro,9o)dro{t) 
(i + diro(t))« 



Hence 



sm,ro,eo) ^ i + d2To{t) ^ i + diTo{t) 

s[l](^;,ro,^o) " ^ ^l + diro(^;) ^^l + diro(^;) 

t'OO 

k{to) < k{oo) = 0(1) / (1 + diwy^dw 

Jo 



Other examples satisfying the integrability condition 2.3 include the inverse Gaus- 
sian core model, and half-symmetric distributions such as the half-normal, half-logistic 
and half-t distribution. 



4 Estimation 



We turn now to estimation of the parameter 6. Let us recall that under the assumption 
of transformation model, the true distribution of (X, 6, Z) is in a class V = {Pe,ri '■ 9 G 
Q,r] & TC}, where 77 represents the triple rj = (F, Gc, A*). 
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To construct an estimator of the parameter 6, we assume that Q is a class of 
probabihty distributions of the variables (X, 6, Z) containing P as a submodel. For 
each {Q,9) e Q x 0, we let Tq q be a monotone function such that 

^Q,0o = To 

if Q = Peo,r)o ^0 — (^0) To, Gc, //) is the true parameter of the transformation model. 

Dropping dependence of this function on the distribution Q, let r„ g be an estimator 
of Tg such that \\Tne — ^e\\oo = oq{1), i.e. the estimate is consistent when observations 
{X, 5, Z) are sampled from a distribution Q & Q. In addition to this we assume the 
following regularity conditions. 

Condition 4.1 Let B{9o,en) denote an open ball of radius £„ and centered at ^o- 

(i) £„ i and ^/rkn T oo- 

(ii) The point Tq = sup{t : EY{t) > 0} is an atom of the marginal distribution of the 
censoring times. 

(iii) The estimate of the transformation satisfies: \/n||rjio— FoUoo = Op(l), limsup„{||r„5) 
ee B{eo,en)} =Op{l) and 

sup{V^||r„, - r„o||oo/[V^I^ -eo\ + i]:ej^eo,9e i?(^o,£n)} = Op{i) 

Examples of estimators satisfying these conditions were given by Cuzick [9], Bog- 
danovicius and Nikulin [5] and Yang and Prentice [23], among others. 

Referring to the notation of section 2, we assume that the function f{u, Z) is of the 
form f{u, Z) — f(r0{u), 9, Z) and make the following regularity conditions. 

Condition 4.2 Let iphe a, constant or a bounded continuous strictly decreasing func- 
tion. For p = 1, 2, 3, let V'p be continuous bounded or strictly increasing functions such 
that < oo and 

poo POO roo 

/ e''^'^\{x)dx < oo, / e''^ip2{x)dx < oo, / e~^'4>^{x)dx < oo . 
Jq Jo Jo 

Suppose that the derivatives of the function i{x, 9, z) satisfy 

\e'ix, 9,z)\< ^Ij{x), \e"{x, 9,z)\< ij{x), \i{x, 9, Z) < ij,{x) 
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The function f{x, 9, Z) is differentiable with respect to x and 

\fix,e,Z)\<M^), \f'{x,9,Z)\<ij2ix). 

We also have 

\gi{x,9,z) - gi{x',9,z)\ < ■max[tps{x),ijj3{x')]\x - x'\, 
\g2ix,9,z)-g2{x,9',z)\<il;s{x)\\9-9'\\, 

where g, = i, f, £" and g^ = 1, f, /', 



The functions ip^ipp^p = 1,2,3 may differ in each inequahty. Note that we do not 
require differentiabihty of / with respect to 9, but only a Lipschitz continuity condition. 

We shall estimate the parameter 9 by solving the score equation [/„(/, 9) — op{n''^/'^), 
where 



-1 " /"TO . 

UnU. ^) = - "E / ^/(^' ^' Zi)M,{du, 9) 



and ^ 

Mi{t, 9) = Ni{t) - I Y,{u)a{V^g{u), 9, Zi)r^e{du) 
Jo 

Here Wf is defined by substituting the estimate Tne into the score function Wf{t, 9, Z) 
of Proposition 2.1. A more explicit form of the score process is given in Section 5. 

Let Eo(/, ^o) = ^o) + Mf. Oo) be given by (2.10) and set 

V{f,9o)= / covo[f-i'^oJ + i'Domu)EN{du), 
Jo 

where ipo = (f* is the solution to the Fredholm equation (2.6). Note that if f{x, 9, Z) — 
£{x,9,Z), then 

V{f,9o) = E^{f,9o) 

+ J covo[f-£''PoJ'](u)(^o + Do[i\(u)EN{du) 
= Ei(/,eo) + S2(/,^o) 



Condition 4.3 (i) The matrix So(/, 9o) is positive definite and the matrix V{f, 9q) 
is non-singular. 
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(ii) The estimate </7„g of the solution to the Fredholm equation (2.6) satisfies 

hmsup„sup{||v9„9||^ : 9 e S(6'o,£„)} = Op{l) and sup{||(^„e - (^olU : e 
B{eo,en)} = op{l). 

The form of the solution (/?o to the equation (2.6) is given in Section 5.3. Therein 
we also verify the condition 4.3 (ii) for the sample counterpart of this equation based 
on an estimator Fno satisfying the conditions 4.1. 

Proposition 4.1 Suppose that the conditions 2.1-2.3 and 4.1-4.3 hold. Then, with 
probability tending to 1, the score equation Un{f,9) = op(n~^/^) has a solution in 
B{Oo,en). In addition, \/n{9 — 9o) converges in distribution to a multivariate normal 
variable Af{0, E(/, ^o)) with covariance matrix E(/, ^o) = (y-^So[y-^]^)(/, ^o)- 

In this proposition the asymptotic covariance matrix of the estimate does not de- 
pend on the estimate of the unknown transformation. In addition, if we choose f — i, 
then proposition 2.1 entails E(/, ^o) — '^o{f: ^o)- 

The second version of this proposition, assumes that a preliminary ^/n- consistent 
estimate 9^^^ of 9o is available. Define 

e^9^^) + vM,o^'^r'Un{f,9^'^) 

Here Vn{f,9^^^) is the plug-in estimate of the matrix Vn{f,9o). Section 5.2 gives the 
explicit form of this matrix. 

Proposition 4.2 Suppose that the conditions 2.1-2.3 and 4.1-4.3 hold. Then ^/n{9 — 
9q) converges in distribution to A/'(0, E(/, ^o))- 



Examples of simple \/n consistent estimators of the parameter 9 were given in 
[5, 6, 7, 9, 12, 23]. 
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5 Proof of Proposition 3.1 



5.1 An auxiliary lemma 

The proof of Proposition 3.1 is based on the following modification of Theorem 2 in 
Bickel et al. [4, p.518]. 



Lemma 5.1 Let B{9o,en) = {9 : \6 — 9o\ < Cn) be a ball of radius £n,£n 0, 
y/nEn — > oo. Suppose that 



(i) v^f/„(^o) ^Ar(0,So(^o)). 

(ii) ^n(^o) ^(^o). 

(iii) The matrices Eo(^o) and ^(^o) are nonsingular. 

(iv) Un{e) - UniOo) = {0- 0ofVM + rem(^), where 



f v^|rem(g) - rem(go)| . ^ A ^ n 

If the assumptions (i)-(iv) are satisfied then with probability tending to 1, the score 
equation Un{0) = op{n~^^^) has a solution 9 in B{9o,Sn) and 

^(9- 9o) ^ 7V(0, [y-^Eo(y^)-^](^o)) 



Proo/ . Let Un{9) = Un{9) - rem(^). We have Un{9o) = Un{9o) because rem(^o) = 0. 
Set 

an = ||/- V^"_^(^o)V;(^o)|| =op(l), 
An = V,-\9o)Un{9o) = Op{n-'/'). 

Finally, define hn{9) ^9- V-\9o)UniO), and put 9l^^ = and 9^^'^ = /inl^n"""^^) for 
m > 1. The condition (iv) implies that for m > 1 we have 

_ ^(0) ^ [J _ (y-iK)(^o)](e-'^ - - ^-^(^0)^7.(^0) (5.1) 
h{et^) - HO'^-'^) = [I- {V-'Vn){9o)]{9^r^ - e-'^) (5.2) 
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Similarly to Bickel et al ([4, p. 518]), (5.1) - (5.2) implies that the mapping /i„ 
is a contraction on the ball = {0 '■ \0 — 6q\ < A„/(l — a„)}. With probability 
tending to 1, C B{9o,en), because An = Op{n^^^'^), an = op(l) and i/n^n T oo. 
It follows that with probability tending to 1, the equation V^^{6o)Un{0) ~ has a 
unique root 9 in Bn satisfying 9 — 9o — V2{9o)~^Un{9o) — Op{n~^^'^). We also have 



Un{9) = U(9) + rem(9) = op(\9 - 9o\ + n''/^) = op(Op(n-V2) + ^-1/2) ^ op(n-V2). 



□ 



5.2 Proof of Proposition 3.1 



Define 



1 



n 



s[l]{t,Tne,9) 



n 



J2y^{tMTne{t),9,Z,), 



i=l 



e[m,rne,9) 



sm,Tne,9) 




Similarly to section 2, we put 



cov[/i,/2](t,r„e,6') 
^r[f]{t,Tne,9) 

^E^^iAr.(i) and set 



ie[fif^]-e{fi]e[gmTne,9), 

6dv[fJ]{t,Tn0,9). 



Let NXt) = n 




p[f, 'Pne](t,Tne,9) 




(5.3) 



cov[/, e']{t, r„,, 9) - var[f ](i, r„,, 9)<pn9{t): 



where (pn9 is an estimate of the solution to the Predholm equation (2.6). 



The score process for estimation of the parameter 9 is given by 
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where 

rt 



Mi{t,e)^Ni(t)- / Y,(u)a{rn0{u),e,Zi)rr,e(du) 
Jo 

and 

Wf{t, e, Zi) = bu{t, r^0, 6) - hiit, r^g, e)ipne{t), 

- [s[l]{t, Fne, e)]-' P Vneit, u)p[f, ip^eKu, r„,, e)NXdu), 

hi{t,Vr,e,e) = /(r,,(0,e,^,)-e[/](i,r,,,e), 
h2i{t, r„,, e) = /(r„,(i), e, z,) - e[e]{t, r„,, e). 

The form of the score process simphfies if we introduce 

"* N,{du) 



Jo 



s[mU,Tng,e)' 



We have 

fTO 



1 " rro 

Un{f,0) = / K(t,r^e,0)-b2i(t,r^e,0)^ne{t)]Ni{dt) 

^ Jo 

[rne-rne]{dt). 



i=l 

TO 







Vne{t, u)p[f, (pn9](u, Tno, 9)NXdu) 

.Jt 

Set M{t) = n-^E2=^Mi{t). Then 



- / Pno - ro](ii)eo[/](ii)ro(dii) + Op(n-^/') . (5.4) 

The score process Un{f,9Q) can be represented as a sum [/„(/, ^o) — Sj=i ^nj(/) ^o)) 
where 



-I "■ /"TO 

Uni{fM = -E/ [&ii(^,ro,^o)-&2i(i,ro,^o)<^o(iMi((^i), 

C/n2(/,^o) = - /"[r„o-r„o]W rVoiu,t)po[f,^o]it)ENXdt) + opin-^/^), 

Jo Ju 

Un3{f,0o) = -/ [rno-ro]{t)po[f,^o]{t)EN{dt)+op{n-'/^), 
Jo 

pro 1 ^ 

UuaU.Oo) = / b„o-¥'o]W-y]&2i(ro(i),^o,i)iV,(cii) = op(n-V2). 

Jo ^ ,_i 
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By central limit theorem, we have y/nUni{f, Oq) ^ -^(0, 9o)). The matrix 9q) 
is defined in Section 2. Further, 



Un2{f,0o) + Un3{f,0o) = - 



■TO 



[r„o - T^o]{du)Vo{u,t) 



Po[f,ipo]{t)EN{dt) 



[r„o - ro](t)po[/, ^omENidt) + op(n-i/2) ^ 
[rno-r„o]((i^) / eo[/](s)ro(ds)Po(s,i) 

J u 



Po[f,^o]it)EN{dt) 



TO 



Uo 

TO 



[fnO-TnoKdu] 



Po[f,^o]{t)EN{dt) 



[r„o - r^o]{s)eo[e']{s)ro{ds)Vo{s,t) po[f,ipo]{t)EN{dt) 
[Tno-Tnompo[f,^omEN{dt) 
- I °[r„o - ro](i)po[/, VomEN{dt) + op(n-V2) . 
Next substitution of (5.3) yields 

f/n2(/,^0) + f/„3(/,^0) = 

- jj ^ (^^ [r„o-ro](^;)eo[f](^;)£;7V(ci^;)^eo[f](s)ro(cis)Po(s,t) 

/ [r„o - ro](s)eo[f](s)ro(ds)Po(s,t)l Po[/,<^o](t)^A^(c^t) 
eo[e]{s)Vo{ds)Vo{s,t) po[f,ipo]it)EN{dt) 
Po[f,^omEN{dt) 



Po[f,ipo]{t)EN{dt) 



+ 
+ 



TO 



TO 



CTO 



dM 



Uo 
^0 ^ dM 



So[l] 

[r„o-ro]Heo[f]^iv(ciM) 



V^o 



^o[l] 



Po[f,<Po]{t)EN{dt) + op{n-'^'') 



Using 



Vo{u,t)-l^ j s^[e]{v)Co{dv)Vo{v,t)^ j eo[e]{v)To{dv)Vo{v,t) (5.5) 

J u J u 
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and Fubini theorem, it is easy to see that the first, the second and the fourth term of 
this expansion sum to 0. The sum of the remaining terms is 



Un2{f,eo) + Usn{f,9o) = - H \ H Vo{u,t)po[f, ^oMENidt) 

Jo so[lJw Uu 



+ op{n-'/'). 



We have ^/n[Un2 + Un3]if,0o) ^ iV(0, E2(/, ^o)), and the matrix S2(/,^o) is defined 
in Section 2. It is also easy to verify that \^[Un2 + t^nsK/, ^o) and ^/nUni{f, Oq) are 
asymptotically uncorrelated. Therefore ^/nUn{f, 9o) => A/'(0, So(/, ^o)), So = + '^2- 

We consider now the expansion of the score process Un{f,0) — Un{f,Oo) for 6 e 
B{eo, En). Set W{e) = f^e - ^ne - Tno + Tno- Then 



w{e){t) = -{e-e') / eo[£]HroW 

Jo 

- [r„e-r„o](t) + Rem(^)(t), 

where the remainder term satisfies 



[Tn0 - rno]{u)eo[i']{u)To{du) 



(5.6) 



sup{ 



^/n\Rem{e){t) - Rem(^o)(t)| . 



■.e^eo,eeB{eo,en),t<ro}^op{i) 



Vn\e -eo\ + i 

For e e 5(6*0, £„), we also have ||r„0 - r„o||oo = op{\e - 6'o| + n-^/^). 
Define 

stiF,(M)/(ro(M), e, Zi)a{ro{u), e^, z,) 



mM,e,Oo) 
e[i']{u,ro,e,eo) 



ns\l]{u,To,9o) 
ns[l]{u,ro,9o) 



ans let 



hnie) 



- J2 / /(ro(w), e, Zi) - /(ro(w), do, z,) 

n Jo L 
-?!/] (u, Fo, e, 0o) + e[/] {u, Fo, ^o)] N^idu) , 

hniO) = -- V r[^'(ro(ti),^,ZO-^'(W,W 

Tl . , .In 



-e[e']{u, Fo, ^, ^o) + e[f ](u, Fq, ^o) Mu)Ni{du). 
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The condition 4.2 implies that Iqn{6), g = 1, 2 is a mean zero square integrable martin- 
gale and Var[^ynIgn{0)] = 0{1){9 — 6*0)^ = = o(l). Using Hocffding projection 
method, we can show that the right hand side of these expressions can be approxi- 
mated by a sum of U- processes of degree k,k < A over Euclidean classes of functions 
with square integrable envelopes Hnk^ k < 4. Under the assumption of the trans- 
formation model, the L2-norm of EHnk is of order 0(£„). Application of maximal 
inequalities for U-processes indexed by Euclidean classes of functions [14] shows that 
Iqn{9) = op{n~^/'^), uniformly in ^ e B{9o,£n)- The details are similar to [12, 13], so 
we omit the proof. 

The score process satisfies Un{f, 9) — Un{f, 9q) — Ej^^/j„(^), where /i„(6') and l2n{9) 
are defined as above and 

"To 



hn{9) = -T r[f(rn0(u),9,Zi)-f{rno{u),9o,Zi) 
-e[f]{u, r„,, 9) + e[/](M, r,o, ^o)]iV,(dM) - 
- E / [^ne - ^no]{u)[f'{ro{u), 9^, Z,) - e[/'](ii. To, 9^)\Nldu) 



n 

i=l 



pro 

(9-9o) / cdv[f,e](u,ro,9o)NXdu) 
Jo 

/•TO 

/ [r„e - Tno]{u)6dv[f,e']{u, To, 9o)N{du) + op{\9 - 9o\ + ||r„e - r„o||oo + n 
Jo 



/•To pro 

= -{9-9o) covo[fJ]{u)N,{du)- [Tne - Tno]covo[f: i']{u)N,{du) 
Jo Jo 

+ Op{\9-9o\+n-'/^) 

1 " /"TO 

hn{9) = — E / [£'(r^6{u),9,Z,)-e\T^o{u),9o,Zi) 

-e[£'](u, Tne, 9) + e[£'](ii, r„o, 9o)\ipo{u)Ni{du) - hn{9) 

^ /°[rne-r,o]H[^"(roH,^o,^«)-e[r](«,ro,M<^oHiv^W 



1=1 



r-To 

+ (^-^o) / 65v[e',e]{u,To,9o)Mu)N{du) 



+ f °[r„e - rno](u)vaj:[l'](u, Tq, 9o)My')N(du) + op(\9 - 9o\ + IKe - r„o||oo + n-^/^) 
Jo 

= {9-9o) covo[e'J]{u)<po{u)NXdu)+ [r„, - r„o]varo[/](K)(^oWA^. W 



+ Op{\9-9o\+n-^/^). 



February 2, 2008 



26 



Combining, 

J2^jn{e) = -{e-Oo) f\oYo[f-i'<Po,i]{u)N{du 

3=1 

/"TO 

- / [Vne-Vn^]{u)po[f,^^]{u)N,{du) + op{\e-eo\+n 
Jo 



The remaining three terms of the expansion are given by 

/"TO pro 

In5{f,e) = - / W{e){du) / Vo{u,t)po[f,ipo]{u)NXdu) 
Jo Ju 

= (^-^o) r eo[t]To{du) ["vo{u,t)po[f,^o]{u)N,{du) 

Jo Ju 

fTo pro 

+ 



pro rro 

/ [r„e-r„o]W / Vo{u,t)po[f,ipo]{t)N{dt) 

Jo Ju 

I °[r„e - Tno]{u)e^y^']{u)V^{du) I ° V^{u, t)po[f, ipo]{u)N{du) 

Jo Ju 



fTo pro 

+ 

+ op(|^-^o| + ||r„e-r„o||) 

/•To 

= -{9 -do) / D4]{u)po[f,^o]{u)NXdu) 



+ 



r[r,, - r,o](«)po[/, ¥^o](/)A^.(c?«) + op{\e - e^] + n-''''), 

Jo 

Ifidif: 0) = - / r„e(?i) , 61) - b2i{u, r„o, 6lo)] [</?„e - ipo\N.{du) 

nj^Jo 

= op(|^^ - ^ol + WTne - Tnoll) = o^{\e - + 

pro 

X I \Vne{u,t)p[f,ipne]{t,rne,0) - Vo{u,t)po[f,^o]{t)]NXdt) 

J u 

[r„o - ^no]{du) 

l-TO 

X / [K(^(M,t)pl/,<^„e](t,r„e,^)-KoKt)pI/,¥^o](t,rno,^o)]iV.(c^^) 

+ rfnO-Vno\{du) 

Jo 

/•TO 

X / [Vno{u, t)p[f, ipno]{t, T^o, ^o) - Vo{u, t)po[f, ipo]{t)]NXdt) 

J u 
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Using (5.4) and (5.5), we find that 

Unif, 9) - UM, Oo) = ^=JUf, 9) = -{9- 9o)VM, 9o) + op{\9 - 9o\ + n-'/^), 

Vn{f,9o) = / coYoif -i'vo,i + i'Domu)N{du). 
Jo 

The matrix Ki(/, 6*0) converges in probabihty to the matrix ^(6*0) defined in Section 4. 
so that Lemma 5.1 completes the proof. □ 

The proof of Proposition 4.2 follows from a similar expansion. The matrix V;, (/. ^^°^) 
can be defined by plugging-in the sample counterpart of the covariance operator in the 
last display. 



5.3 Verification of the condition 4.3 

We have shown in [12] that the equation (2.10) simplifies if we multiply both sides of it 
by 7^(0, Let V^(0 =mO"'V'(i),5[/](i) =7^(0, i)-iL'[/](i) and p[/,-L'[/]](i) = 
ViO,t)p[f,-D[f]]{t) Set 

c(t)^ [ V{0,u)-^dC{u), b{t)^ [ V{0,ufB{du). 
Jo Jo 

Multiplication of (2.7) by V{Q,t)-^ yields 

i/j{t)+ k{t,u)ij{u)h{du)= j k{t,u)p[f,-D[f]]{u)EN{du) , (5.7) 



where the kernel k is given by k{t, u) = c{t A u). The square integrability condition 2.3 
is equivalent to the assumption that k{to) — c{u)b{du) is finite. The solution to the 
equation (5.6) is given by 

ip{t)= rA{t,u)p[f,-D[f]{u)EN{du), (5.8) 
^0 

where A(t, u) = A{t, u, —1), and A(t, u, A) is the resolvent corresponding to the kernel 
k. The solution to the equation (2.6) is given by 

cp{t) = -D[f]{t) + r A{t,n)p[f, -D[f]]{u)EN{du)V{0,u)V{0,t) (5.9) 
^0 
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Prom [12], A{u,t) is given by 

^^"'^^^ ' ^'-''^ 

where for s <t, the interval functions "i/olsjt) and ^'i(s,t) are defined as solutions to 
the Volterra equations 



^^0(5,^) = 1+/ c{{s,u])b{du)'ifo{u,t) ^1+ '$o{s,u-)c{du)b{[u,t]) 

J{s,t] J{s,t] 

= c{(s,t])+ c((s,u])b{du)^i(u,t) ^c{(s,t])+ ^i(s,u)b(du)c((u,t]). 

J(s,i\ J(s,t] 



Define also 



^'2(5,^) = 1+1 b{[s,u))c{du)'^2{u,t) = 1+ I '^2{s,u+)b{du)c{{u,t)) 

J[s,t) J[s,t) 



*3(s,t) = b{[s,t))+ I b{[s,u))c{du)^3{u,t) = b{[s,t))+ ^3{s,u)c{du)b{[u,t)) 

[s,t) J[s,t) 



Then 



*o(s,t) = 1+/ '^i{s,u)b{du) = 1+ c{du)'^3{u,t+) , 

J{sA J{s,t] 

*i(s,t) = / ^Q{s,u-)c{du) = I c{du)-^2{u,t+) , 

J(s,t] J(s,t] 



^2{s,t) = 1+/ b{du)^i{u,t-) = 1+ ^3{s,u)c{du) , 

Jls,t) Jls,t) 

*3(s,t) = / ^2{s,u)b{du) = / b{du)'^o{u,t-) . 

J \s,t) J \s,t) 



If To is an atom of the survival function P(X > t), then = 0,1,2,3 form 

bounded monotone increasing interval functions. In particular, \E'o(s,t) < exp k,{tq) 
and 'i/i{s,t) < 'i/o{s,t)[c{t) — c(s)]. If Tq is a continuity point of the survival function 
P{X > t) and k(to) < 00, then ^o{s, t) < exp k,{to) for any < s < t < Tq, while the 
remaining functions are locally bounded [12]. 



Next suppose that the transformation model holds with {0,r) = (^OiTo)- We 
assume that the estimate Tno satisfies the conditions 4.1 and show the natural plug-in 
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estimator ip^e of (5.9) satisfies the conditions 4.3. Define 

Cne{t) = I Vne{O,u)~'^Cn0{du) 
Jo 

hne{t) = [ Vne{0,uf^r[i']{u,Tne,9)NXdu), 
Jo 

where Cne and Vne are defined as in Section 5.2. Then the sample analogue of the 
equation (5.9) reduces to a system of linear equations which can be solved by inverting 
a bandsymmetric tri diagonal matrix [11, 12]. 

Denote by ^n9,j the sample counterparts of the interval functions ^j, j — 1, 2, 3, 4. 
Using Fubini theorem we can show that 

[*ne,o - ^o](s,^) = / '^n9,o{^,u-)cne{du)[bne-b]{[u,t]) 

J{s,t] 



+ / [cne-c)]{{s,u])b{du)^o{u,t) 

l{s,t] 



/ '^ri0,o{s,Ui-)CneidUi)[bn0 " 6] ([^i, W2))c(rf«2)*3(^i2, *+) 

J S<tll<«2<t 

+ / '^ri0iis,ui)bneidui)[cne - c]{{ui,U2\)b{du2)'ifo{u2,t) 

J S<Ul<U2<t 



and 



[*nM-*l](«'^) = [Cne - C]{{s,t]) 

+ / '^ne,l{s,u)bn0{du)[Cn0 - C]{{u,t]) 
J(s,t] 



'{s,t] 

+ / [^nefi-'^oKs^u-Xdu) 

J{s,t] 



Under assumptions of the condition 4.1, we have c„ = sup{|c„£i — c\(t) : 9 G 
B{eo,e„),t < To} and b^ = sup{|6„e - b\{t±) : 9 e B{9Q,en),t < Tq} -^p 0. 
In addition, for = 0, 1, we have 

limsupsup{v['„gg(0,r) : 9 G B{9o,en)} 

n 

< [c(ro) - c(0)]'^exp[/^(To)](l + 0^(1)) = 0^(1) 

Hence 

f = sup{|^„e,o-^o|(s,t) :0<s<i<T,^eS(^o,£n)} 
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< 26Ae,o(0, To)^^o(0, To) + 2cAe,i(0, r)^'3(0, r) 
sup{|^^„0,i - ^i\is,t) ■.0<s<t<r,9 e B{9o,en)} 

< 2c„*„e,o(0, To) + rc(To) 

and both terms converge in probability to 0. 

Let A„e(s, t) be defined similarly to (5.10) The preceding calculations, imply that 
sup{A„0(s,t) : s,t e [0,To],9 e B{9o,en)} = Op(l) and sup{|A„e - A|(s,t) : s,t e 
[0,Tq],9 G B{9o,Sn)} = op(l). Verification that the sample analogue (fine of the equa- 
tion (5.9) satisfies the conditions 4.3 can be completed using Gronwall's inequalities 
given in [12] and integration by parts. 
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